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We study a class of strongly interacting many-body fermionic systems in 2 + 1-dimensional non- 
relativistic conformal field theory via holography. The 5-dimensional charged black hole with asymp- 
totic Schrodinger isometry in the bulk gravity side introduces parameters of background density and 
finite particle number into the boundary field theory. We propose the holographic dictionary, and 
realize a quantum phase transition of this fermionic liquid with fixed particle number by tuning 
the background density /3 at zero temperature. On the larger j3 side, we find the signal of a sharp 
quasiparticle pole on the spectral function A{k,uj), indicating a well-defined Fermi surface. On the 
smaller /3 side, we find only a hump with no sharp peak for A{k,uj), indicating the disappearance of 
Fermi surface. The dynamical exponent z of quasiparticle dispersion goes from being Fermi-liquid- 
like 2; ~ 1 scaling at larger /3 to a non-Fermi-liquid scaling 2 ~ 3/2 at smaller /3. By comparing 
the structure of Green's function with Landau Fermi liquid theory and Senthil's scaling ansatz, we 
further investigate the behavior of this quantum phase transition. 
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I. 



INTRODUCTION 



Experimental physics on strongly interacting non-relativistic many-body bosonic and fermionic systems develops 
rapidly on account of the progress of controUing ultra cold atoms [1]. The dynamical exponent z of those 
microscopic states is usually z ^ 1, along with other well-known condensed matter systems such as non-Fermi 
Uquids metals from heavy fermions and high Tc superconductor(see [2] and reference therein), beyond the paradigm 
of Landau Fermi hquid theory. In this strongly coupling regime, the traditional perturbative field theory study 
has been challenged, while holography, specifically Anti de Sitter space and conformal field theory (AdS/CFT) 
correspondence[3], becomes a powerful alternative. Holographic methods have shown some success in the study 
of certain strongly interacting fermion systems [4-7], with the emergence of Fermi surface and non-Fermi liquid 
behavior. (See [8-10] for reviews on the holography applied to condensed matter physics.) However, thcre are at 
least two major shortcomings to bridge this success to ultra cold atomic systems. On one hand, the dual field 
theories of these study are asymptotic conformal and relativistic, which are quite different from the non-relativistic 
nature of ultra cold atoms. On the other hand, this AdS/CFT setting bears no parameters identifying the particle 
number, mass or density Spectrum, because of these parameters are not good quantum numbers in the relativistic 
theory. Tuning physical parameters such as particle number or doping density becomes important in the absolute 
zero temperature, where purely quantum fluctuations can drive phase transitions, known as quantum phase 
transitions[12]. 

Substantial works in the literature had contributed to understand non-relativistic conformal field 
theory(NRCFT)([13-15] and reference therein), as a renormalization group(RG) fixed point of the non- 
relativistic systems. And its gravity dual theory had been proposed, with solutions of zero temperature pure 
Schr6dinger(Schr) geometry[16, 17], finite temperature black holes[18-20], and charged black holes[21, 22]. There 
has been some pursuits on studying fermions in this asymptotic Schrodinger geometry [23-25]. However, to our 
knowledge, the holographic study of Fermi surface from strongly interacting fermions with NRCFT background has 
not yet been explored in the literature. Our paper is aimed to bridge this gap and tackle the two aforementioned 
shortcomings. 

Schrodinger black hole in the bulk gravity theory, on one hand, realizes an asymptotic NRCFT background with 
the dynamical exponent z = 2 for the boundary field theory naturally. On the other hand, Schrodinger black hole 
provides the particle number(or the mass operator[26]) M from the gauge invariant ^-momentum: M = £ — qA^\g 
[26] and background density /3 to the non-relativistic boundary field theory. Our approach is similar to AdS/CFT 
set- up [5], considering a Dirac fermion field in the probe limit under a charged black hole spacetime, the Green's 
function can be read from the asymptotic behavior of Dirac fermion field in the UV of the bulk side. Here 
we also propose the holographic dictionary of real-time retarded Green's function for fermions in Schr/NRCFT 
correspondence[27], analog to the work of [28] for AdS/CFT case. For convenience, we name these classes of 
strongly interacting non-relativistic fermionic liquids under asymptotic NRCFT background as Schrodinger Fermi 
liquids. 

The paper is organized as follows. Firstly we discuss the charged Schrodinger black hole solution and its Dirac 
fermion equation of motion, to introduce our notations in Sec.II. We then provide our holographic dictionary [27] 
analog to the setting of [5, 28] in Sec.III. In Sec.IV, we demonstrate the zero temperature nearly ground state of 
this fermion system shows a sharp quasiparticle peak in the spectral function - the evidence of Fermi surface, with a 
non-Fermi liquid dispersion relation. We compare our spectral functions A{k,uj) of Schrodinger Fermi liquids with 
Landau Fermi liquid theory and Senthil's scaling ansatz[29]. In Sec.V, we show the evidence of a quantum phase 
transition, by tuning the background density /? but fixing particle number at zero temperature. On the larger /3 
side, we find a well-defined Fermi surface. On the smaller f3 side, we find only a hump with no sharp peak for 
A(k,uj), indicating the disappearance of Fermi surface. The dynamical exponent z of the quasiparticle dispersion 
goes from Fermi-liquid-like scaling z = 1 at larger /3 to larger z(~ 3/2) non-Fermi liquid at smaller /3. Finally, we 
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conclude with some remarks and open questions in Sec.VI. 



II. SET-UP: DIRAC FERMION FIELD IN A CHARGED SCHRODINGER BLACK HOLE 



Based on the holography, a quantum field theory of finite charge density can be mapped to a charged black hole 
of a gravity theory [5]. The U(l) charge of Schrodinger black hole induces finite charge density to the boundary 
field theory, meanwhile breaks the non-relativistic conformal invariance. Thus, we only have 'asymptotic' NRCFT. 
Bcfore discuss the details of bulk gravity theory, it will bc hclpful to introduce generic labels for a large class 
of NRCFT(with charge and mass densities) we will study. We characterize our 'asymptotic' NRCFT by five 
parameters, (A, M, fiq, P, T)[26]. Two parameters, the conformal dimension A and the mass operator M = £—qA^\g 
from the gauge invariant ^-momentum, define the boundary NRCFT in a universal sector. The U(l) charge chemical 
potential /zq and other relevant terms from the current J^(such as charge density pQ and mass density pm) in 
NRCFT is mapped to U(l) gauge field of the bulk gravity. Background density /3 is introduced by Schrodinger 
black hole through NuU Melvin Twist(or TsT transformation) [18-22]. The physical way to interpret this /3 could be 
the density of doping background, or an analog of interaction strength t/U of Hubbard model^. Temperature T of 
the boundary theory is given by the black hole temperature Tbh- Notably, the conformal dimension of NRCFT here 
depends on the mass operator M, which is quite different from CFT. More peculiarly, the conformal dimension for 
spinors has an extra rnzL^ split, as already been noticed in [23, 24]. In the foUowing we denote dimension d as the 
spatial dimension of xi^X2, ■ ■ ■ ,Xd- Thus the bulk spacetime of Schrodinger asymptotic as Schr^+s (distinguished 
from the AdS(i+2), where as the corresponding boundary theory of Schr^+a is NRCFT^+i (CFT^+i for AdSd+2)- 
We summarize the conformal dimensions of spin-0 boson[26] and spin-1/2 fermion operator([23, 24]) in the foUowing 
table. 



asymptotics 


AdSd+2 


Schrd+3 


scalar conformal dimension 






spinor conformal dimension 




A± = ^ ± ^'(m ± 1)2 + M2 



TABLE I: conformal dimensions of CFT and NRCFT for spin-0 boson and spin-1/2 fermion 



A. Charged Schrodinger black hole 



We focus on c? = 2, 5-dimensional(5D) Schrodinger black hole Schrs in the bulk and 2+lD NRCFT3 on the 
boundary. Let us briefly go through our set-up for the charged Schr5. In string frame, the metric is, 

d4tr = ^ {~fdr' + dy' - p\'f{dr + dyf) + ^(dx? + dxj) + (ILl) 

where R is curvature radius. By converting to light-cone like coordinates, t — f3{T + y), £, = + y)i ^nd 

switching to Einstein frame for the later use of holographic dictionary. 



rf4» - K-y^ l^"^ [C-^ ~ r'f)dt' + - f)de + (1 + f)dtd^) + "^.{dxl + dxl) + 
where f{r) ^ 1 + ^ - (r^ + ^ j ^,K(r) = i+B^rHi-f(r)) - ^^'^ charged black hole supports a gauge field. 



r4'-"-V'y— l-|-/32r2(l-/(r)) • 

A^Adr, A.-^(l-^), (II.2) 



^ Indeed /3 has dimension [length]^, so the "background density over area" with correct dimension should be defined as /3 ^ Wg simply 
name /3 as background density for convenience. 
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Adopted in [26] notation for later use, A = ^dt — jSArd^ = Atdt + A^d£^, with 

At = ^lQ+ PQ , A^=Mo + PM r-^ (11.3) 

By holography[26], fig is identified as U(l) charge chemical potential, pq,pm are the charge density and mass 
density, Mo is related to the mass operator hy M = £ — qA^\g = £ — qMo with £ as the ^-momentum. The 
temperature of the black hole is given by identifying the inverse of the near-horizon Euclidean periodicity of 
boundary time coordinate i, 

r.»-r?kfi-SV (II.4) 



Our interest of study is the boundary field theory at zero temperature, which corresponds to the extremal black 
hole with Q — •\/2rg. This being said, all the numerical analysis contained in this paper pertains to zero temperature 
only. At zero temperature, the charged black hole Schrs horizons degenerate, meanwhile the near horizon geometry 

becomes AdSa x , with^ ds"^ = -e^'^^'V^Ads^ + ^AdS2(*Ve^) + r^dx^/R\as2' ^^^h i?AdS2 = i^- 

B. Dirac fermion 

To probe the fermionic response of the boundary theory via holography, we proceed to solve the Dirac fermion 
equation in the bulk curved spacetime of the charged Schrodinger black hole. The action is 



Sotrac = y ^ x^J-gEinii'{ei'T'"Di_, - m)^), (II. 5) 

and its equation of niotion(EOM) is {e^V^V^ — m)ip = O, with covariant derivative 

= 9,, + iryaec.f^[^^ r^] - iqA^, (II.6) 

where gamma matrix of flat tangent space {r'',r''} — 2rj°-^ , vielbeins relates flat tangent space to curved 
spacetime, ff^iye-^eg" — 77.^. The spin connection is w?^ — e^d^j^ey + V^^e^e^ , F"^^^ are the Christoffcl symbols. 

We choose the specific vielbein'^, 

* 2r^fB ' t 4r2/B ' ? ^ ^' ^ 2/32(/-l) ' ^ ' 

= Xh/r, eg' = i^h/r, e/ = K'^r^f. (II.8) 

with other unwritten components of are zeros. To simplify the Dirac equation calculation, here i? is a function of 
r chosen ensuring the coefhcient of r*r«r^ in the EOM is zero. The boundary behavior of B{r — > 00) is a constant 
B, times l/r. We choose B^ to be 1. Asymptotically, B(r -> 00) ~ i + |^ + '^'-^£5"''''' + -i6+60/3^^+24ff^+27ff« ^ 

-^^ 256r9 + ■ ■ • . The near-horizon behavior of B is Bh/{r — 1), where Bh is a constant. Given 

= 1, we can numerically solve this equation to find Bh- In 5D spacetime, each of matrices has 4x4 
components, we choose to express them as follows: 

r, = 



^ see Appendix A for details 

^ For convenience, we rescale the coordinates to set R = tq = 1 from now on. 
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(:/»)• --(ro)' -Ho") 

(Ji are Pauli sigma matrices, / is identity matrix. Yt — /3(r^ + rj^), = (— + r.y)/(2/3). Rewrite the 4-component 
Dirac spinor field '\\) as: 

where ø+ and are two-component spinors. This (—5.9'^'^)^^^* factor eliminates a term in the Dirac equation, 
which is simplified to 

where u and v are hnear combinations of the vielbein components e-*,e-^, e~ and e/: 

" = + 9^*) - A + ~ gA^) f/3e/ + ^^ef] ) (IFU) 



v^K-^Å{Lo-r qAt) [pe} - ^e|) + (£ - qA^) + ^e|) 



(11.12) 



By rotational symmetry of the boundary theory, we wiU work on the case fci = O and set k2 — k from here on. 
We will write (f> — (ø+,ø_)"^, also its ø+ = (j/+, z+)'^ and ø_ — (y_, z^)"^ in the component form. 

III. GREEN'S FUNCTION FROM HOLOGRAPHY 

A. Holographic dictionary 

We study the fermionic response of the boundary theory, by probing the Dirac fermion field in the bulk 
spacetime of Schrodinger black hole. The holographic dictionary of source-response relation can be set up by 
reading the boundary action of Eq.(II.5). From [28, 31], the variation of bulk action induces a boundary term 
Sg = Jgji^ (Pxd^^y—gg^^'ipip. Thcrefore the relation between bulk field and its conjugate momentum are, 

n+ = n_ = y=^^+ (iiFi) 

We can identify the source(x) and response(C') from boundary(or UV) behavior of bulk field(V'±) and 
momentum(n±), from the holographic dictionary, 

eM-Sgrav[^,^]{r ^ CX))] = (exp[ J d'^+'xixO + dx)])QFT (111.2) 

The source x ^-nd bulk field are related by, 

X = lim rT^-''±'(/' (IIF3) 

r— f 00 

the response O and momentum Il± are related by, 

O = - lim r^^-^Il (III.4) 



r—yoo 



where h'± = y {m ± ^)'^ + {i — qMoY generically^ , analogue to the result of [28]. Here we show only the standard 
quantization(corresponding to source A and response D of [28]). The alternate quantization(corresponding to 



^ specifically equal to (m it 1/2)2 _|_ _|_ qQfj")2 jjj q^j. charged Schrodinger black hole case 
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source C and response B of [28]) can be done in the same manner^. The Green's function Gr is related to the 
ratio of O and x- 

We now study the Dirac equation (11.10) in boundary UV asymptotic hmit to extract O and x from the coefficients 
of lp and n, or equivalently related to (/)+ and ø_ at r — >■ oo. In this hmit, (11.10) becomes 



r^/fdr T mK W'± + 7^ P± + -P± H Pt ± — <^2 ± — cri U>t = O 

V /\2p r r r r J 

where 



C = §^ (-4'ZQ(1 + /?') + 5(£ + gQ/3)/33 - (£ + gQ/?)g'/3' - 8/?u;) , P± = 

(/)+ = Si r'^+^s (Al + A2r-2) r-'^+^s (ai +a2r-2) + 83 r''- + 5 (Bi +B2r~2) +R2 r-"- + 5 (/3i +/32r-"') + • ■ • 

(III.7) 

= Si r''+ + 5 (Cl + C2r-2) + ^-"^+ + 5 (71 +72r-2) ^ 33 r"— 5 (25, + D2r-^) + R2 r-"^— 5 (<5i +(52r-2) + . . . 

(III.8) 

here v± = ^(m ± 1/2)2 _^ (^ + gQ/3)2. Each of Si, S2, Ri, R2 is a r-independent one-component multipher, as the 
coefficient of the spinor®. Thcrc is a projection relation between the two-component spinors^: 

SiCi = — (^±^P_SiAi, S2i?i^— (^±^P,S2i.i, (III.9) 

2l3{v+^m+\) 2/3(iy_-m+i) 

R-iTi = 7^, — — TT^-R-i"i' ^2l3i = — — ^P+R25i, (III.IO) 

2/?(-i^+ + m + i) 2/3(-i^_ - m + i) 

We now apply our holographic dictionary to identify the source and response from (III. 7), (III. 8). To read the 
boundary value, in the foUowing we take r — >■ 00 as UV limit. Consider the leading behavior of contribution, 

^- = i-ggn-'^U- ^ r-^_ c, s^Cr-^^-'/^ (III.ll) 

which corresponds to the source 

X- = lim r5-''+V^_ ~ SiCi (111.12) 
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For the alternate quantization, 

X = lim r^-^^ip (III.5) 



r— f 00 



the response O and momentum Il± are related by, 



O = lim r^i-T^n (III.6) 

r — ^oc 



® When doing numerics for this field redefinition, it is important to keep subleading term C2 in the Si r'^+ ^2(... + C2»" ^ + ...) series, 
since this C2 term dominates Ri r~'^+~ 2 ai when uj^ > ^, which is indeed our case in the numerical study. Thus here we keep the 
expansion for all four sets of solutions to the subleading orders. 

Each of Al, j42, Cl , C2, Ol, «2i 7l! 72> Bl , _B2, i^i , D2, /32i <5i , ^2 is a two-component spinor. List above totally there are sixteen 
two-component spinors. The spinors Ci, C2, 71 and 72 are in the nuU space of P4., the spinors Bi, B2, /3i and /32 are in the nuU space 
of P_. There are four independent sets of bases, each basis as a solution of Dirac EOM; the first set contains Ai, A2, Ci, C2 and its 
subleadings, the second set contains 01,02,71,72 and its subleadings, the third set contains B\, B2, D\, D2 and its subloadings, the 
fourth set contains /3i , /32 , <5i , (52 and its subleadings. 
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is proportional to Si. The momentum field n_ is 
which corresponds to the responsa O-, 

0_ = lim r''+~5n_ ~ Ri«! 
We take the asymptotic constant term on the UV boundary. is proportional to Ri . 



On the other hand, we can go through the same logic again, though consider the leading behavior of > 
bution, 



which corresponds to the source x+j 

X+ 

X+ is proportional to 82- The momentum ficld n_|_ is 

n+ = -V^^- = -i-ggn'^U- ^ 



lim r2 

r^oo 



^4 



which corresponds to the response 0_|_, 



O, 



■ lim 



in. 



R2(Si 



(III.13) 



(III.14) 



contri- 



(III.15) 



(III.16) 



(III.17) 



(III.18) 



0+ is proportional to R2. Now we derive 81,82 are identified as sources, Ri,R2 are identified as responses for 
this standard quantization. A similar argument works for the alternative quantization by taking the subleading 
terms of the bulk field and its conj ugate momentum, we leave this detail in the Appendix C. 

In addition to our above holographic dictionary, we provide another intuitive argument on identifying source and 
response. We notice the boundary action V'^A, due to the F^. form, it couples the first componcnt of the spinor to 
the third component of "ø, while couples the second component of to the fourth component of "ø- Both two-point 
function or {pip shows r^'^+ or r'^'^- forms[23], which pairing r'^+ to r~'^+ , or r'^' to r~'^- . Thus it is consistent with 
our earlier dictionary set-up. We will perform a more constructive comparison with [23] and a pure Schrodinger 
Green's function computation via our dictionary in Appendix D. Lastly, we are aware that the detailed construction 
of Schr/NRCFT holographic dictionary involves with nontrivial holographic renormalization[23, 24, 32, 33]. Our 
works here only present the spirit of source-response holographic dictionary following Ref.[5, 28]. Apparently, 
the rigorous holographic regularization and renormalization for spinors are the future next step to justify our 
Identification of Green's function. 



B. Source and response from UV expansion 



To extract the data of source and response, we define a converting matrix Cv as a function of r(see Appendix 
B), and a set of functions 8i(r), 82(f), Ri(r), R2(r) satisfies 



^) - (ø+(r) 0-(r))^ = (y+{r) z+{r) y^(r) z_(r))^ - Cv • [Si(r) Ri(r) S^{r) B.^{r) 



(III.19) 



By this field definition, neatly Si(r), S2(r), Ri(r), R2(r) approach to 81, 82, Ri, R2 as r 00. Due to 
projection, we find the spinors have the properties ci+ = C2+ = 71+ = 72+ = 61- = ^2- = = P2- = 0. To 
deal with the standard quantization, from the lesson of Sec.IIIA, we choose ci_ = 1 to compute the first set of 
bases, ai_ = 1 to compute the second set of bases, 61+ = 1 to compute the third set of bases, (5i+ = 1 to compute 
the fourth set of bases. Each of four independent bases in (III. 7) and (III.8)(equivalently in B.2, see Footnote.7), 
can be determined by a free parameter, thus totally four free parameters. Now the four free parameters for four 
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independent bases are 81,82, Ri, R2- Argue from Sec.IIIA, for the standard quantization, the source terms are 
8iCi,S2-Bi, with their corresponding responsa terms Riq;i,R2^i respectively. Our choice of spinor = (0,1) 
and its coupled spinor aj — (0:1+ , 1) justifics that coefficient 81 is exactly a source and Ri is its response. Our 
choice of spinor Bf — (1,0) and its coupled spinor — (l,(5i_) justifies that coefficient 82 is exactly a source 
and R2 is its response^. 



C. IR behavior and the In-falling boundary condition 

To determine the near horizon initial condition of Dirac equation, here we deal with IR behavior and solve the 
in-falling boundary condition at zero temperature, Q = y/2. Consider the small e expansion of the equations, where 
r = 1 + e. The equation for B{r) becomes ^ = —l/e- Thus, we take the behavior of B{r) near horizon as -B/i/e, 
where Bh is another constant. limr_>.i(r — l)B{r) = Bh For the later use, we define 



2/? 



-/3a; 



(III.20) 



Cj is the coefRcient of r in the exponent dependence of Eq.(II.9). We find that the behavior of /, u,v in the IR is 
given by 



/->12e2, u 



-[/, V 



-V 



(III.21) 



4^3 V 2V^Bh 



with U 

wave function ansatz are 



and V 



4^3 



+ 



Dirac equation near horizon and its infalling 



2V3 



Cd{U±Va3 



(111.22) 

t>± CX e' (III.23) 

The exponent of wave function (j)± is chosen to be + sign, in order to combined with Eq.(II.9) to be e~*'^'^+'"/(i^'^) 
infalling into the black hole^. The infalling condition is obtained by plugging Eq.(III.23) into Eq.(III.22), where 
the subscript H stands for values at the horizon. 



b+\H = {U + Va3)cj)^\H 



(111.24) 



The infalling condition for spinors has two linear independent choices, the first set is ø-.i = (1,0) thus = 
{U + V,0), and the second set is (j)-^2 = (0, 1) thus (j)^_2 = iO,U — V). Therefore this gives two independent sets 
of initial conditions at horizon for Si(r), Ri(r), 82(f), R2(?'), which we introduce one more upperindices 1,2 to 
distinguish the first and the second sets: 



[Slir)Illir)Sl{r)Rlir)f\H ^ G^~'\h 



[Slir) Klir) Slir) Rl(r)f \„ 
More conveniently in matrix form. 



8(r) 



81 (r) Slir) 
Slir) Slir) 



H 



H 



C. 



R(r) 



Rl(r) Rlir) 
Ri(r) Ri(r) 



^\h ■ (u + V O 1 O 



O u-v o 1 



For the alternate quantization, we should alternatively choose ai_ = 1 to compute the first set of bases, 7i_ = 1 to compute the 
second set of bases, cii+ = 1 to compute the third set of bases, = 1 to compute the fourth set of bases. From Sec.III A, the source 
terms are 82^1, SiAi; with their corresponding response terms R2/3i, Ri 71 respectively. Our choice of spinor Af = {ai_|_, 1) and its 
coupled spinor 7^ = (0, 1) justifies that coefficient Si is exactly a source and Ri is its response. Our choice of spinor Df = (1, 
and its coupled spinor = (1,0) justifies that coefficient S2 is exactly a source and R2 is its response. 
The 12 factor appears here origins from the near horizon geometry AdS2. 
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D. Green's function 

Green's function of the boundary theory is defined to be the ratio between source matrix S(r) and response 
matrix R(r). Thus we define G(r) based on R(r) = G(r)S(r), and evaluate G(r) at r — > oo, to read the 2x2 
matrix Green's funetion G of the boundary theory^", 

G = Hm G(r) hm R(r)S(r)-i (111.25) 

We derive the EOM of G(r) in the bulk gravity(see Appendix B), and solve this EOM with the initial 
condition: G(r)|/f = R(r)|^f • S(r)^^|H- to obtain physical results of Eq.(III.25). 

IV. SPECTRAL FUNCTION 

We now equip with the holography tool developed in Sec.II,III. The original questions driving our interests 
are: what is the nearly ground state of this fermionic system under asymptotic NRCFT background at zero- 
temperature? Will there be a Fermi surface? Will Fermi surface coUapse, destabilized by tuning non-tcmperature 
parameters (such as background density /?)? Will this realize certain quantum phase transition of fermionic 
liquids? With the holographic dictionary for Green's function, we proceed to study these questions. We focus on 
as zero temperature phase. 

The spectral function A(k,uj) of this boundary system can be determined by the imaginary part of Green's 
function. In 2 + 1 D boundary theory with 2x2 matrix G, we should take eigenvalues of G, namely, 
Im[G(r oo) eigenvaiuesl ■ Similar to [5], among two eigenvalues, one of eigenvalues Im[Gi] is picked for de- 
tailed studies in our analysis, the other eigenvalue is less relevant to the physics for experimental comparison, such 
as ARPES data. Below we will simply abbreviate Im[Gi],Re[Gi] as ImG,ReG. 

A. Fermi surface 

This non-relativistic fermionic system has five parameters, conformal dimensions ^±{i^±)-, temeprature T, chem- 
ical potential(of background) fig, particle number eigenvalue or mass M, and background density /3. We first study 
the background density a,t (3 — l/\/2 at T = 0. The gauge-invariant mass operator M = £ — qM^ — £ + qQ(5 is 
fixed to be 1/10. The Dirac fermion charge q = 1, its mass is chosen to be to = 1/10, nonzero value in order to 
avoid scaling dimension i'± degeneracy and extra logarithmic term in UV expansion. We find there is a sharp pole 
on Im[G'] at ujp = 0.8984, fc^^^ = 1.3169, indicating a stable quasiparticle like excitation at Fermi-momentum kp. 
This indicates a well-defined Fermi surface. Normally the location of Fermi surface on oj-axis is shifted by chemical 
potential /j,, one redefines ui = lj — fi thus uj = O has the Fermi surface. In our case, is shifted by the presence 
of ^-momentum £, this can be realized from the faet that the location of Fermi surface is determined mainly from 
the low energy IR physics. which in the bulk gravity corresponds to near horizon region. Thus, instead of using 
boundary time coordinate t and its coupled conjugate energy cj, we identify the near-horizon time coordinate r and 
its conjugate energy uj. When ui — {jg + f3uj)\i^p = O, namely ujp = —£/{2f3'^), its value indicates the pole location 



There is no extra factor multiplied with G(r) for this Green's function, becausc in our dictionary sources and responses are related 
to the coefficients of two-component spinors, instead of spinors itself. 
In principle, 

A(k,ui) = --lm[G] (IV.l) 

TT 

in our Schrodinger holography case there is a subtlety. When we relate the d + 3-D gravity theory to a d + 1-D boundary theory via 
Eq.(III.2), the constant factor J = needs to be absorbed into J d'*+^a;(xC' + Ox), this gives an extra constant factor for source 
field or response field. Namely, the A{k,(jj) may be different from Im[G] with another extra constant factor. This factor should be 
important when justifying spectral density sum rule, ^ I^^c ^ik,^)duj = 1. The exact value of our Im[G] is less informative, only 
the relative height of Im[G] has physical indication. In addition, note somc proposes to study the trace part[6], Im[Tr[G(r — > oo)]]. 
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of a Fermi surface. Denote k± = \k — kp\, we find near the quasiparticle like peak has scalings, 



Im[G(a;*(fc_L), kj_)] 



Z Cl 
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FIG. 1: At /3 = l\/2, (a)Imaginary part of Green's function, with a pole. (b)Real part of Green's function switches sign 
Crossing zero valne at a specific uj near the peak of Im[G']. 
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FIG. 2: (a)At /3 = l/\/2, the imaginary part of Green's function, Im[G] as a function of ui and k. A sharp quasiparticle- 
hke pole at ujf = 0.8984, = 1.3169 indicates a well-defined Fermi surface. The pole indicates infinite lifetime stable 
quasiparticle at fc^. Notice the main branch of dispersion goes into ui < ujf and k > kp, & hole-like excitation. While in [5], 
their main branches of dispersion goes into cj > uup and k > kp, a. particle-like excitation. (b) The dispersion relation traces 
the pole into four branches on {k,u}) plane in a large scale. The asymmetric behavior indicates a particle-hole asymmetry. 



In Fig. 1, wc show the imaginary part Ini[Gi] and real part Re[Gi] of Green's function, see the location of Fermi 
surface indicates a pole in Im[Gi] and switches the sign of Re[Gi]. The main branch of dispersion goes into lo < ujp 
and k > kp, which is a hole-like excitation The result is difFerent from [5], where their main branches goes into a 
particle-like excitation with uj > ijjp and k > kp. In Fig. 2, the dispersion relation shows particle-hole asymmetry 
in large scale, though close to {ujp,kp), it gives a unique dynamical exponent z. 
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B. Comparison to Landau Fer mi liquid theory and Senthil's scaling ansatz 



To better understand the physics of Green's function G{k,uj), we now study the functional form of G{k,uj) in 
terms of two different classes. Both classes hold under general arguments. The first class is Landau Fermi liquid 
theory, which holds for weak coupling system, where the free fixed point is still a good description of the system. 
The second class is even more general based on scaling ansatz for non-Fermi liquid theory and critical Fermi surface, 
proposed by Senthil[29]. In Landau Fermi liquid(LFL) theory, the retarded Green's function is of the form, 

11 Z 

^^^''^^ " c^-a-S(cj,Å:) " cj-(^fe+ReS(fc,cj))-zImE(fc,w) ~ (uj - LUp) - (k + ^ ^"^^'^^ 

S(a;, k) is the particle irreducible retarded self-energy. = e/j — fi, is the excitation around the original chemical 
potential. The condition + Re Wi?) = ujp to define renormalized Fermi-momentum kp. The final form 

is obtained by expanding ^j. + ReS](fc,a;) around {kp^uip)^ with the definition of quasiparticle residue Z, with 
Z^^ = 1 — ReS(fc = kp,uj — iop), also = (fc — kp)Z -^{S^kf + KeT,(kp,ujp)) = v{k — kp), and quasiparticle 
decay rate I/t^ = — 2ZImS(/c,CLi). The specific LFL form we use to fit our Green's function is 

^ = n ^ f , ^^ ■ f . (IV.5) 

[—{lu — LUp) — V{k — kp)) — tJ[LL!) 

with our quasiparticle self-energy ansatz as 7(0;) = k{uj — cjp)", where k is some real constant, LFL has n = 2. We 
will take general n for fitting ansatz. We flip the sign of (w — tup) to have a hole-like dominant excitation as Fig.2 
suggests. 



The scaling ansatz proposed by Senthil[29] based on general arguments, has the form at T = O, 

G = co{k-kp)-'^F„Cjl^^) (IV.6) 

for better fitting we will be forced to choose cq and ci their values on two sides ui > uip and lo < ujp differently to 
reflect particle-hole asymmetry. 



In the following subsections, we present our Green's function data for k < kp, k = kp and k ~ kp, and fit these 
data by LFL and Sentil's ansatz. The main messages of our fitting(Fig. 4, Fig. 6, Fig. 8) of LFL and Sentil's ansatz 
to our data(Fig. 3, Fig. 5, Fig. 7) are: 

(1) Sentil's ansatz generally has better agreement than LFL fitting for our IniG data. 

(2) Our Re G data are sandwiched by the LFL fitting with LFL fitting with n = 2 and Marginal Fermi liquid (MFL) 
with n = 1[11, 30], which likely implies that our Schrodinger Fermi liquids can be a closer description between LFL 
and MFL theory with 1 < n < 2. From our quasiparticle self-energy ansatz as 7(0;) = k{lo — ujp)'^ and quasiparticle 
decay rate l/r^ 7(w), this may suggest Schrodinger Fermi liquids has shorter life time and larger decay rate 
l/^fe ~ (w — wf)^"^ than LFL l/rk ~ (w — LOpY close to Fermi energy LUp. Compared to LFL, the quasiparticle 
description of Schrodinger Fermi liquids is less robust. 

(3) Near the pole location, we have not found a promising fitting for Sentil's ansatz for both sides oi lo > ujp and 

LO < LOp. 
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1. k < kp 




FIG. 3: fc = 621/500 < kp. Blue curve is for ImG, red curve is for Re G. 



For k < kp, our scaling ansatz is 

cojk - kp ) 

Cl (k—kf) 



. , - (IV.; 
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(a) <^ (b) <^ 

FIG. 4: fc = 621/500 < kp, both (a)(b) with three fitting curves: (l)LFL Eq.(IV.5) with n = 2 in gray-dashed, (2)LFL 
Eq.(IV.5) with n = 1 in purple-dashed, (3)scaling ansatz form Eq.(IV.8) in black-dotted-dashed. Co,Ci,7o are chosen to be 
positive but their valnes on two sides u > cup and uj < ujp are chosen differently. 



For ui < uip under k < kp, thc term inside logarithmic becomes negative, which we choose the complex logarithm as followini 

co(fc - kp) 

(uj—ujp) 



(IV.7) 
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2. k> kf 



k>kf 




FIG. 5: k — 8/5 > kp- Blue curve is for ImG, red curve is for Re G. 
For k > kp, our scaling ansatz is 

Co(fc-fcf)"" 

iog( ,;(rC)^j -»7o 
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FIG. 6; fc = 8/5 > kp, both (a)(b) with three fitting curves: (l)LFL Eq.(IV.5) with n = 2 in gray-dashed, (2)LFL Eq.(IV.5) 
with n = 1 in purple-dashed, (3)scaling ansatz form Eq.(IV.lO) in black-dotted-dashed. co,ci,7o are chosen to be positive 
but their valnes on two sides lo > lof and o; < cjf are chosen differently. 



For LO > LOp under k > kp. thc term inside logarithmic becomes negative, which we choose the complex logaritlim as foUowini 

co(fc - kp) 

(uj—ujp) 



, , (IV.9) 
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3. k = kf 



k=kF 



ImGi ,ReGi 




FIG. 7: k = kp. Blue curve is for ImG, red curve is for Re G. 




FIG. 8: fc = kp, both (a)(b) with two fitting curves: (l)LFL Eq.(IV.5) with n = 2 in gray-dashed, (2)LFL Eq.(IV.5) with 
n = 1 in purple-dashed. At fc = kp, scaling ansatz of both forms Eq.(IV.8) and Eq.(IV.lO) runs into trouble. We have not 
yet found a good fitting. 



LFL fitting must be symmetric respect to w = ivp at k = kp for ImG, however the ImG from Schrodinger Fermi 
liquids is not symmetric along w — ujp. This is the major difference. 
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V. FERMIONIC QUANTUM PHASE TRANSITION 

Schrodinger black hole introduce two additional parameters to the non-relativistic conformal background, other 
than the parameters occurred already in asymptotic AdS spacetime[5]. The first parameter is the mass operator 
M(particle number eigenvalue), the second one is background density (3^^. In the zero temperature phase, a natural 
question arises: what happened to the Fermi surface if background density /3 is tuned? We should fix M while 
varying /3 In Sec.VA, we encounter again the phase with Fermi surfaces, graduahy tune down /3 near 1/2, we 
find a critical point(or critical line) in Sec.VB. Smaller /? shows that Fermi surface collapses and then disappears 
in Sec.VC. Altogether may indicate a quantum phase transition of fermionic liquids. 

A. Well-defined Fermi surface {/S > f3*) 

When /3 > 1/2, an obvious peak appears in IniG, seeFig.9. Analytically the peak should approach 5(fc — fci?,CLi — 
LUp) at (kpjUJF)- However, the numerical value lm.G{kp,ujp) cannot really be infinite. What we find is that the 
peak lmG{kp,ujp) values in this region /? > 1/2 depends on the finite IR cutoff. The smaller the initial cutofF 
e — r — I, the larger the numerics lTaG{kp,ujp) at the peak grows. This is a sign for the supposc-to-bc infinite 
pole. The pole of ImG indicates a well-defined Fermi surface. At larger f3, the pole and nearby region on (/c, w) 
plane develops much sharper. Notice that the pole shifts to larger ujp and smaller kp hy decreasing /?. 

B. Near the quantum critical point (/? ~ 

As (3 approaches in the range between 1/V2 and 1/2, we find the lmG{kp,ujp) peak becomes insensitive to IR 
finite cutofF e, the peak values are lower for smaller f3, shown in Fig. 10. The stable peak value indicates there is 
no (5-function like pole on (fc,w) plane. By tuning /3 to smaller value, the Fermi surface gradually coUapse and 
disappear. We interpret the physics as: 

(3 > 1/2, liRG{kp,ojp) ~ Z5{k - kp, uj - ujp) + finite terms, where Z 7^ 0. (V.l) 
/3~l/2, lmG{kp, LOp) Z5{k - kp, uj - ujp) +fLmte terms, with Z ^ O (V.2) 

Since the Z goes to zero at finite /3 ~ 1/2, we suspect it is not a smooth crossover behavior. We expect a quantum 
critical point (3* (or quantum critical line) slightly larger than /3 = 1/2, and smaller than /3 = 1 /V^. For convenience, 
we denote /3* ~ 1/2, as in Fig. 10(b). We do not numerically determine /?* due to the computational limitation. A 
more detailed scan near the peak at /3 ~ 1/2 may determine the exact value of (3*. 

C. Fermi surface collapse and disappearance (/3 < /3*) 

At smaller /3 < 1/2, we find no sharp peak but only a smoother hump. Unlike Mott insulator for Mott transition 
[29], we do not have gaps opened up in A(k, uj). Indeed A{k, uj) does not dip to zero in this phase. There is no non- 
analiticity in A(k,uj) to pin point kp. This shows it is still a gapless phase but with Fermi surface disappearance. 



Usually a quantum phase transition is tuned by a dimensionless coupling[12], in our case we can define the coupUng = fj^/iQ 
with fj,q fixed in our case. Since the Hamiltonian description of boundary theory is unknown, we schematically tune /3 to address 
the same physics as tuning g^. 

Here we tune /3 with various values 1/16, 1/4, 1/2, l/\/2, 1, 2, 8. The other parameters should be fixed. We choose M = £ + 
fixed to be 1/10, T = 0{Q = y/2), /^±{v±) is fixed by m = 1/10 and M. Among all the five parameters of the system, the remained 
parameters is subtle, which is iiq = Q/(2/3). In our numerics, we choose to fix g = 1, in this case it seems like chemical potential 
liq varies while /3 is tuned. One may argue that a resolution is considering = q^iq where fiq still allowed to be fixed while q 
compensates to be adjusted correspondingly. This resolution seems to fix the (chemical potential) energy to add a fermion of charge 
q into the system. However, we should aware that in any case the chemical potential /jg is indeed varied. In addition, the 'real' 
chemical potential to set the scale of Fermi energy iip is not merely as in [5] only alone. In our Schrodinger system, the Fermi 
energy should be identified by the coordinate u), the Fermi energy is set by oj + = (/3(tj + qAt) + {1 — 9^{)/ (2/3)) |a = Co- This 
shows that = q(^(5At — Aj/(2/3))|g = O is independent of /3. Therefore, Fermi energy /ip = O is already fixed. We choose to fix 
the charge q of fermionic contents, instead of varying fermion charge g to fix the energy fj^q of inserting a fermion 
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We show a series of ImG plots by varying /3 in Fig.9,10,11. Note that the vertical axes for ImG shows no tick 
marks, we only show a landscape scanning through many slices of ImG. Each slice of ImG(/c,a;) has a fixed w, 
and scanning k valnes. Each slice of ImG has been shifted vertically for a clear vision of the landscape. As in 
Footnote.ll we had discussed the size of the compact circle modifies ImG to the physical value of A{k,uj). 
Thcrefore, here the exact valuc of ImG is immaterial, only the relative height of ImG matters. 
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FIG. 9: ImG of Schrodinger Fermi liquids for /?>/?*: (a)/3 = 2 (b)/3 = 1 
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FIG. 10: ImG of Schrodinger Fermi liquids: (a)/3 = 1/V2 (b)/3 = 1/2 ~ /3*. Near P* . 
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FIG. 11: ImG of Schrodinger Fermi hquids for ^ < ^*: (a)/? = 1/4 (b)/3 = 1/16 
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D. Evolutions of dynamical exponent z and Fermi- momentum kp under tuning background density /3 

Here we study the evolutions of the dynamical exponent z, Fermi energy ujp^ Fermi-momcntum kp while tuning 
/3 in Fig. 12, Fig. 13. When P P* , there is no good quasiparticle description for the system, so in which case 
(wf, kp) means the (w, k) coordinates of the highest peak in spectral function, z just means the dispersion reading 
from the branches structure around the highest peak. 




4 5 8 0.0 0.5 1.0 1.5 2.0 2.5 

P (b) kAP) 



FIG. 12: (a)The relation between z and P for six data points are: /3 = 1/4,2 = 1.593, /3 = 1/2, z = 1.299, /3 = l/\/2, z = 
1.143, /3 = 1,2 = 1.050, ^ = 2,2 = 1.021, /3 = 8,2 = 1.009. (b)For various /3 = 1/4, 1/2, 1/V2, 1, 2, the Fermi-momentum 
and energy {kp, cup) are (0.000030303,2.020), (0.8762,1.212), (1.317,0.8984), (1.789,0.656), (2.8304,0.3266) 

In all of the data above, a ~ 1.00, thus our data foUows the general relation z > a and 2; > 1 as Senthil's 
argument[29]. As f3 increases, z goes close to 1. Tentatively it suggest a more Landau Fermi liquids like behavior 
at large /? limit. Though from the spectral density fitting, we find the imaginary part of quasiparticle does not 
obey 7(w) cx and Im[G(fc,w)] near kp is not symmetric respect to ujp. These two features are distinct from 
LFL. The large f3 limit is at most a close cousin of LFL. 

In Sec(IVA), we showed oop = —t/{20^). In the case of tuning /3 while fixing the mass operator I — qMo = 
I + qQP = M, we expect ujp{(3) = — (M — qQf3) /{2/3'^). We show this power law fitting agrees with our data in 
Fig. 13(a). On the other hand, the Fermi-momentum kp requires better understanding of UV physics, we do not 
have a fitting here. 



2.5 FT 
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FIG. 13: (a)The location of the Im[G] peak at ljf varies respect to j3. The numerical data points compared to the dashed 
curve fitting (jJf{P) = —(Af — gQ/3)/(2/3^). (b)The location of the Fermi-momentum kp varies respect to /3. 
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VI. CONCLUSION AND OPEN QUESTIONS 

In summary, we have studied a class of strongiy interacting non-relativistic fermions under asymptotic NRCFT 
background in 2 + 1 D. We make some efForts to deal with the aforementioned two shortcomings of AdS space. 
Firstly, our model has a better realization of non-relativistic properties of many body systems, and we have 
observed the well-defined Fermi surface of Schrodinger Fermi liquids. Secondly, by tuning the background density 
j3 with fixed particle number M, we realized a fermionic quantum phase transition as Fig. 14, where on larger /3 
side shows a sharp Fermi surface, while on smaller /? side shows Fermi surface disappearance. We find Senthil's 
scaling ansatz generally has better fitting agreement than Landau Fermi liquid with our non-Fermi liquids. Based 
on quasiparticle self-energy scaling, we argue the quasiparticle description of Schrodinger Fermi liquids has shorter 
lifetime and is less stable comparing to LFL. 

We leave some questions for future directions: 

(1) Quasi-particle residue Z may be regarded as the order parameter for the quantum phase transition. How 
does Z in Eq.(V.l) behave near quantum critical point (or line), what is the order of phase transition? We have 
not answered these questions due to numerical power constraint. It will be illuminating to understand whether 
Schrodinger Fermi liquids shows discontinuous Ist, or continuos 2nd order or higher order transition, and the 
possibility to realize similar phase transitions proposcd in [29]. 

(2) Notably the charge U(l) symmetry and particle number U(l) symmetry are unbroken in our probe limit. 
Fermi surface and gauge-gravity duality relation are mentioned in [34], especially the relation between a global 
U(l) symmetry and the existence of Fermi surface. How does our system realize a quantum phase transition 
with Fermi surface disappearance without breaking global U(l) symmetry or translational symmetry? Bring our 
attention is an earlier work[26], we consider a toy model of bosonic system under asymptotic NRCFT, where 
a U(l) symmetry is broken by condensed boson helds around Schrodinger black hole. Bosonic quantum phase 
transition is likely found there at low temperature phase as Fig. 15. It is unavoidable to ask whether these two 
phase transitions in Fig. 14, Fig. 15 have any similar nature. On one side, /3 > /3* of Fig. 14 shows a conducting 
phase with Fermi surface with unbroken U(l) symmetry; il > 17* side of Fig. 15 shows a metallic state with 
unbroken U(l) symmetry. On the other side, /3 < /3* of Fig. 14 shows Fermi surface disappearance; $7 < il* side 
of Fig. 15 shows a superfluid state with broken U(l) symmetry. Though the two systems have similar asymptotic 
NRCFT background, beware the two systems are rather different. The bulk side of fermionic model has a charged 
Schrodinger black hole, where the gauge field is chosen to be fixed, the Dirac fermion is in a probe limit. On 
the other hand, the bulk side of bosonic model in [26] has a neutral Schrodinger black hole, where both the 
gauge field and bosons are in a probe limit. And the notions of these two parameters /?, il can be quite different, too. 

(3) It will be illuminating to address more about the phase with Fermi surface disappearance. There is no 
superconducting gap opened up in the spectral function of the phase without Fermi surface. Specifically we do not 
introduce any pairing term(such as Yukawa spinor-scalar pairing) in the bulk action, so it is not a superconducting 
phase^^. We only can suspect that tuning j3 from large to small effectively implies tuning fermion interaction from 
weaker coupling to stronger coupling - from more Fermi-liquid-like(z ~ 1) to non-Fermi liquids(z > 1) to strongiy 
correlation smear the A{k,uj) discontinuity into continuity near cup, kp. Whether one can understand more about 
the nature of this Fermi surface disappearance, we leave this for future study. 

(4) It is of considerable interest to perform rigorous holographic regularization and renormalization for spinors to 
justify holographic dictionary of Green's function, foUowing [23, 24, 32, 33]. Indeed, we find a subtle issue already 
that the conformal dimensions A of Schrodinger spinors have h'± with peculiar m± 1/2 dependence, distinct from 
AdS case[5, 28]. In AdS case, degrees of freedom rely on two sectors of two component spinors(in [5, 28] notation, 
D and A for the standard quantization, B and C for the alternate quantization) . In Schrodinger case, degrees of 
freedom rely on two sectors of z/+(as Si, Ri) and two sectors of i^_(as S2,R2) shown in Eq.III.7,III.8, but with 
only an one-component(the coefficient of the projected spinor) degree of freedom for each sector. One may wonder 
whether it is possible to maintain two-component spinors structure as AdS case to construct the Green's function in 
Schrodinger case, however we find that this approach will sacrifice the distinction between the standard quantization 



A set-up along [35, 36] is the next-stcp toward boson-fermion interaction, with a superconducting state under NRCFT background. 
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and the alternate quantization[37], which is unreasonable. A more rigorous justification is desirable for future study. 

(5) Our model here is for a 2 + ID fermion system. It will be important to study a system in 3 + ID which 
may exhibit fermions at unitarity[15, 16]. It will also be interesting to explore dynamical exponcnts other than 
z = 2. Indeed gravity duals of finite density systems with asymptotic Schrodinger isometry for d 7^ 2, z ^ 2 are 
known[19, 38, 39]. 

(6) The Schrodinger black hole at zero temperature has finite entropy, which implies that our theory do not 
describe a unique ground state but an ensemble of low energy states. Moreover, it has been pointed out that 
discrete lightcone quantization and (3 deformation from the parental AdS black hole[18-20, 40, 41] causes peculiar 
free energy scaling F ^ —T'^/jj? for the system. It will be interesting to know whether bosons or fermions with a 
full consideration of spacetime back reaction can change the physics of our study, especially the IR AdS2 geometry. 
The hope is that our study still characterizes the physics of nearly ground state nature of asymptotic NRCFT 
system with fermions. 



S 

No FS FS 



FIG. 14: The zero T quantum phase transition diagram of Schrodinger Fermi liquids. A phase with Fermi surface(FS) in 
P > 13* . A quantum critical region near 13 'r^ j3* (it is undetermined yet whether (3* is a critical point or critical line, and 
unknown whether the transition is 1-st order or higher order). A phase without Fermi surface in /3 < /3*. Note to tune a 
dimensionless couphng gp, we can define gp = PJJIq with /iq fixed. 




Superfluid Metal 

FIG. 15: Bosonic quantum phase transition is likely found at low temperature phase in [26]. On large background density (SI > 
Q,*) side, the phase is in metallic state with unbroken U(l) symmetry. On smaller background density(f2 < f^*) side, there 
shows a superfluid phase with broken U(l) symmetry. Note to tune a dimensionless coupling gn, we can define gn = ^/ y/JJKi 
with fixed. 
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Appendix A: AdS2 scaling 



The near horizon geometry of charged Schrodinger black hole can be obtain by taking r = ro + e, with the 
coordinate redefinition: 



-1/3 



1 + P^rl 



2/), V 
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I2 ' 



(A.l) 
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K^ = {l + P^rl)-\ 

With RAdS2 = ^ — R ^ VrT^' ^^^"^^ ^^^^2 xM-^. The gauge field near horizon is: A — Ardr = 
{ArKQ^'^)dT + {—Ar(3^rQKQ^^\/l2)dy, with Ar ~ 7^~3 £ = — One can solve the Dirac equation, with this 

AdS2 background and supporting gauge field. We take the AdS2 rescahng as in [22], send t ^ t/X and e — > Ae 
with A — > 0. In this case, the Dirac equation near the AdS2 boundary becomes. 



(-^r,9, + wj^Tr + CfTr + CyTy + C,-,r,, - m)V' = o (A.3) 

JiAdS2 ^^AdS2 



Rewrite the above in terms of two sets of two-component spinors and "square" the operator to make it a second 
order differential equation, we find the "AdS2" scaling dimension is the exponent v oi tp (x e~^^'^. 



V = RAds, Jrn^ + C?~ Cy2 - C,~^2 (A.4) 



with Cf = -Z g,gg^^^ , Cy = zi?AdS2 VT2i^o + 2^), = iRAAS,k^/V2Kl'^ . 



Appendix B: Numerical Set-Up for Green's function 

For the foreseeing convenience, we define a matrix Dm, which satisfies Dirac equation Eq.(II.lO), (f>' = dr(t> — Dm 



Dm = 



1 



O 



y ikjr —u — V 



—u — V —ik/r \ 

—ik/r —u + v 

-mK-^/^ O 

O -mK-^l^j 



(B.l) 



We define a converting matrix Cv as a function of r as 



C. 



„1/2 



(ci_ + C2-r-2) (7i_ + 72-7-2) 



(di_+d2-r-2) 



(,5i_+<52-r-2)^ 



(B.2) 



and a set of functions Si(r), S2(r), Ri(r). R2(r) can be defined from Eq.(III.19). This field-redifinition Si(r), 
S2(r), Ri(r), R2(?') goes to 81,82, Ri, R2 at r -J> 00. 

The EOM of G(r) in the bulk gravity is G'(r) = R'(r)S(r)-i-G(r)S(r)'S(r)-i. Apply Eq.(B.l) and Eq.(III.19), 
then S'(r) and R'(r') can be simplified in terms of linear combination of 8(r) and R(7') 



8i'(r) Ri'(r) 82'(r) Ra'w] = (Cv~V(r))' = (Cv~' Dm Cv - Cv-'Cv') [8i(r) Ri(r) 82(0 R2(? 
The matrix Dg{r) ee Cv^^ Dm Cv - Cv^^Cv' simplifies EOM to. 



(B.3) 



G'(r) 



"^92,1 


^92,3' 


+ 


"^92,2 


Dg2A 








_^S4,2 


Dg4,4 



G(r) - G(r) • ( 



Dgl,l 


Dgl,3 


+ 


'Dg2,l 


Dg2,3 


_Dg3,l 


Dg3,3 




_Dg4A 


Dg4,3_ 



G(r)) (B.4) 



Numerically we solve this bulk EOM Eq.(B.4) of Green's function with the initial condition to obtain physical 
results of Eq.(III.25). 
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Appendix C: Holographic dictionary for the alternative quantization 

We also walk through the similar game for the alternative quantization. Consider the subleading behavior of (p- 
contribute as a source field, then 

V'- = i-gg^-'^U- ^ c SaCir--^/^ (Cl) 

which corresponds to the source 

X- = hm rS-'^-T/'- ^ S2D1 (C. 2) 

r— >-oo 

X_ is proportional to S2. The momentum field II_ is 

II- = V^gg^^lJ+ = {-ggn'^U+ ^ r^cl^+ ~ S2Bir''-+5/2 + R2P,r-'-+^^^ + ... (C.3) 
which corresponds to the response 

a_ = hm r^—^U^ ~ R2/3i (C.4) 

r— >C30 

O- is proportional to R2. On the other hand, we can go through the same logic again, though consider the 
subleading behavior of ø+ contribute as a source field, then 

= i~99n-'^U+ ^ r-^<t>+ ^ Si^ir^^+'S/^ (C.5) 

which corresponds to the source x+: 

X+ = lim ri-''+^+ ~ SMi (C. 6) 



X+ is proportional to Si. The momentum field n_|_ is 

n+ = -V^W^iJ- = -i-grf'^^- ^ -r'P- ^ -SiCir''++5/2 _ R^^^^--++5/2 + . . . (C.7) 
which corresponds to the response 0+, 

0+ = - lim r'^+-in+ ~ Ri7i (C. 8) 

r— ^00 

0+ is proportional to Ri. Now we again derive Si, S2 are identified as sources, Ri,R2 are identified as responses 
for this alternative quantization. 



Appendix D: Pure Schrodinger Green's function at zero T zero density 
1. two-point correlators 

Here we compare the pure Schrodinger two-points function at zero T zero density of Ref . [23] with our formula- 
tion-'^^. Specifically, in d = 2, they show, 

(V'M(a:,t)V^A/(0,0)) (xr2''+ (D.l) 

Here we crosscheck our analysis indeed matches theirs. The boundary action is dtd^cp'x^—gg^'^'^t^, with 
V~99^^ =J^j by plugging our UV boundary expansion in Eq.(III.7), Eq.(III.8), with V' = {—g9^^)~^^'^4> — ^~^</') 
we arrive '4''4'\dM ^ r~^(j)(j). Notice Tr in cfxf) coupling Ist to 3rd component of (j), meanwhile coupling 2nd to 4th 
component of ø. With leading piece in ø is the 4th component of ø, which is Si r'^+^s . This couples to the leading 
order of 2nd component of 4>, which is Si r'^+~2a2. Neglect other factors and coefficients, 

JdM 

FoUow the same logic along [23], one can obtain the same scaling and same form of two-point correlator 

{'^M{x,t)i^M{OM- 



In [23], their e coordinates are inverse of our r, also their d = 3 is our d = 2 case 
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2. Pure Schrodinger Green's function from response over source 



Solve the Dirac's equation in zero T zero density Schrodinger spacetime ds^ — —r'^dt'^ + 2r'^dtd£,+r'^dx^ + dr'^/r'^, 
thc bulk fieids have the following form: 



i'-= /+(k,r)r^v^++ 



g+(k,r)r^C/++ r-^i^_,^(k/r)T/- + 



-X,_(k/r)C/+ + 



/_(k,r)r^v^- 



.g_(k,r)r4C/- (D.3) 
^if_^_(k/r)C/" (D.4) 



where d is the spatial dimension of the quantity k = \J + fc^ is a coordinate-invariant form of momentum. 
We denote U'^ ^V^ as two components spinor bases for v± series(the first two columns), and denote C/^,V^^ as 
two components spinor bases for —v± series(the last two columns), K^ji^jr) is modified Bessel function, solution 
oir'^dlK^C^Ir) + rdrK^{\Llr) - {iy^jrf ^v'^)K^i^lr) =0. The expansion of if^^ (k/r), /±(k, r), 5±(k, r), V+, V-- 
near r — > oo are 



K,^ (k/r) = 2-1+"^* (k/r)-'^±r(j.±)(l + »(r'^)) 
/±(k,0 = + 0(r-2)) 



±j/+ + rn - 



5±(k,r) = -i- 



^^(±"-)fcT-r-^±-(l + 0(r-^)) 



(D.5) 

(D.6) 
(D.7) 



V'+ = Ar''+-— (l + 0(r-2)) + Sy.'^-- — (1 + 0(^.-2)) + 



V'- 



^(1 + 0(r-2)) + /Sr-'^— ^(1 + 0(r-2)) 

(D.8) 

Cr"+-'^ (1 + 0(r-2)) + Dr''—^{\ + 0(r-2)) + 7r-''+-^ (1 + 0(r-2)) + (^r-"--"^ (1 + 0(r-2)) 

(D.9) 



Notice "0+ = {~99rr) ^4'+ ^-nd ■!/;_ = {~99rr) '^4>~i can compare this expansion respect to Eq.(III.7), 



Eq.(III.8). The expansion matches, with the projection constrains on the spinors: 



C 



2/3(z^+ + m+ i) 



-P^A. B 



2/3(t._ -TO + i)^+^' 



7 - 



Y-P-a, P 



2/3(-i^_ -m + i) 



(D.IO) 
(D.ll) 



In the case of Charged Schrodinger black hole for Eq.(III.7), Eq.(III.8), the projection relation Eq.(D.lO), 
Eq.(D.ll)'s € is replaced by ^ + qQI3 for Eq.(III.9), Eq.(III.lO). 
There are extra constraints on two-component spinors , : 
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2£ 



or equivalently, 



2£{kt'T^) 



2£(fc'^r^) 



k2 



(D.12) 



(D.13) 



where fc^F^ = iTt — ujT^ + kx^x- By identifying source and response based on our holographic dictionary, we 
have the response and source matrix, 



R 



R2 R| 

Si SI 

§2 S2 



2-''+-ir(-j/+)fc'^+(v-)^^'^^' 
2-''--^r{-i'^)k-'- {u-)i^^^^ 



2-''+^'^T{-v+)k''^ 
2-''--^V{-v^)k''- 



Hik'^V^T(^V-)'i^^'> 

K^fc^r^r^c/-)^^) 



^k--{T^U+)\ 



(1,2) 



(D.14) 



(D.15) 
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Green's functioii G = RS has the similar structure as [23]. 
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